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Abstract 

The purpose of the present effort is threefold. Firstly, it is shown that 

there exists a principle, that we call Kinetical Interaction Principle (KIP), 

underlying the non linear kinetics in particle systems, independently on the 

picture (Kramers, Boltzmann) used to describe their time evolution. Secondly, 

the KIP imposes the form of the generalized entropy associated to the system 

and permits to obtain the particle statistical distribution, both as stationary 

solution of the non linear evolution equation and as the state which maximizes 

the generalized entropy. Thirdly, the KIP allows, on one hand, to treat all the 

classical or quantum statistical distributions already known in the literature 

in a unifying scheme and, on the other hand, suggests how we can introduce 

naturally new distributions. Finally, as a working example of the approach 

to the non linear kinetics here presented, a new non extensive statistics is 

constructed and studied starting from a one-parameter deformation of the 

exponential function holding the relation f{—x)f{x) = 1. 
PACS number (s): 05.10.Gg, 05.20.-y 
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I. INTRODUCTION 



In the last few decades there has been an intensive discussion on non conventional clas- 
sical or quantum statistics. Up to now several entropies with the ensuing statistics have 
been considered. For instance, in classical statistics, beside the additive Boltzmann-Gibbs- 
Shannon entropy which leads to the standard Maxwell-Boltzman statistics, one can find in 
the literature, the entropies and/or related statistics introduced by Druyvenstein |]T]J^, Renyi 
§, Sharma-Mittal §, Tsalhs §, Abe §, Papa §, Borges-Roditi §, Landsberg-Vedral 0, 
Anteneodo-Plastino Frank-Daffertshofer |TT], among others. On the other hand, in the 
literature we can find, besides the standard Bose-Einstein and Fermi-Dirac quantum statis- 
tics and/or entropies, the ones introduced by Gentile ]12|], Green [|l3l, Greenberg-Mohapatra 
1^ , Biedenharn [|l^, Haldane-Wu |ll6| , p^ , Acharya-Narayana Swamy ||T8|, Buyukkilic- 



Dimirhan [|T9l etc. This plethora of entropies poses naturally some questions. 

A first question is if it is possible and how to treat the above entropies in the frame of 
a unifying context and from a more general prospective, in such a way to distinguish the 
common properties of the entropies, from the ones depending on the particular form of the 
single entropy. 

A second question is if it is possible to obtain the stationary statistical distribution of 
the various non linear systems in the frame of a time dependent scheme. Fermion and boson 
kinetics were introduced in 1935 by Uehling and Uhlenbeck |2^. On the other hand, the 
non linear kinetics associated with the anomalous diffusion has been considered in the last 
few decades in several papers by the mathematicians and by the physicists of the condensed 
matter. After 1995, in the frame of the Fokker-Planck picture, the anomalous diffusion has 



been linked with the time dependent Tsallis statistical distribution [p| -|30| and the kinetics 



of the particles obeying the Haldane statistics and the quon statistics |^ has been 
considered. 

The problem of the non linear kinetics from a more general point of view has been 
considered only in 1994. In ref. [p3| it has been proposed an evolution equation (eq. (7) 
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of the reference) describing a generic non linear kinetics. Subsequently some properties of 



this kinetics in the frame of the Fokker-Planck picture has been studied in ref. [31| and in 



the frame of the Boltzmann picture in ref. [Q. Finally, the kinetics described by non linear 



Fokker-Planck equations has been reconsidered recently in ref. |Tl| , |35[] . 

It is well known that the formalism used to describe the time evolution of a statistical 
system, depends on the picture used to describe the system. For instance, for a particle 
system interacting with a bath, we can study its time evolution in the phase space in the 
frame of the Kramers picture (Fokker-Planck picture in the velocity space). Besides, for an 
isolated system, we can study its time evolution in the phase space adopting the Boltzmann 
picture. 

A third question which arises at this point is if the entropy of a system, or its stationary 
statistical distribution, depends and how on the particular picture used to describe the 
system. 

A fourth and last question is if it exists a principle underlying the time evolution of 
the system, in the two pictures. Obviously if this principle exists, it must define both the 
entropy and the stationary statistical distribution of the system. 

The present paper is concerned with the above questions. Its principal goal is to show, 
that exists a principle in the following called Kinetical Interaction Principle (KIP), which 
governs the particle kinetics and imposes the form of the entropy of the system independently 
on the particular picture used to describe the system. Within the two pictures and in 
a unifying context, the H-theorem is proved and the form of the the generalized entropy 
together with the stationary statistical distribution for a generic non linear system are 
obtained. 

The paper is organized as it follows. In Sect. II, we introduce the KIP underlying the 
kinetics of a particle system, without regard if it interacts with its environment or if it is an 
isolated system . In Sect. Ill, we study the nonlinear kinetics in the Kramers picture implied 
by KIP. In particular, after writing the evolution equation of the system we obtain its entropy 
and also its statistical distribution both as the stationary state of the evolution equation and 



by using the maximum entropy principle. The stabihty of this equihbrium distribution is 
also studied. In Sect. IV, we examine an isolated particle system and describe its nonlinear 
kinetics governed by the KIP in the Boltzmann picture. In particular the equilibrium and 
the stability of the system are studied. In Sect. V, we consider in an unifying context some 
examples of already known classical and quantum statistical distributions, in order to test 
and highlight the utility of the formalism here developed. In Sects. VI and VII, we consider, 
just as a working example, a new statistical distribution, the K-deformed distribution, which 
arises naturally in the frame of this formalism. This distribution is obtained according to 
the KIP both as stationary solution of a nonlinear evolution equation and by using the 
maximum entropy principle. In Sect. VIII, we consider some concrete physical systems 
where the K-deformed distribution can be adopted. Finally in Sect. IX, some concluding 
remarks are reported. 

II. A PRINCIPLE UNDERLYING THE KINETICS 

Isolated systems: Let us consider an isolated system composed by identical particles. 
We make the hypothesis that the system is a low density gas so that we can describe it 
by a one particle distribution function. The interaction of a particle in the site r = {x,v) 
where the particle density is / = f{t,r), with a second particle in the site Vi = {xi,vi) 
where the particle density is /i = f(t,ri), changes the states of both the particles. After 
the interaction we find the first particle to the site r' = {x',v'), where the particle density 
is /' = fit, r') and the second particle to the site r'^ = {x[, v[), where the particle density is 
f[ = f{t, r[). We postulate that the transition probability from the state where the particles 
occupy the sites r and ri, after interaction, to the state where the particles occupy the sites 
r' and r'^ is given by 

^(^t, r^r', r,^r[) = Tit, r, t\ n, r[) i{fJ')i{hJ[) . (1) 

The first factor in (|I|) is the transition rate which depends only on the nature of the two 
body particle interaction. Then this factor is proportional to the cross section of the two 
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body interaction and doesn't depend on the particle population of the four sites. Being the 
system composed by identical particles, the second and third factors are given by the same 
function and differ only on the arguments. The factor 7(/, /') in (|I|) is an arbitrary function 
of the particle populations of the starting and of arrival sites. Eq.(0) takes into account 
two body interactions and for this reason the function 7(/, /') must satisfy the condition 
7(0, /') = because, if the starting site is empty, the transition probability is equal to zero. 
The dependence of the function 7(/, /') on the particle population /' of the arrival site plays 
a very important role in the particle kinetics because can stimulate or inhibite the particle 
transition r — > r' in such a way that interactions originated from collective effects can be 
taken into account. The condition 7(/, 0) 7^ requires that in the case the arrival site is 
empty the transition probability must depend only on the population of the starting site. 
We note that for the standard linear kinetics the relation 7(/, /') = / holds. 

Systems interacting with a hath: We consider now the case of a particle system interacting 
with its environment which we consider as a bath. The particle which transits from the site 
r to the site r' is now the test particle while the one that transits from the site Vi to the 
site r'^ is a particle of the bath and the factor entering into the transition probability is 
indicated with 7 and depends on the nature of the bath. If we address our attention 

to the test particle and after posing 

iy(t,r,r')= / rf2Vid2v;T(t,r,r>i,r;)7(/i,/0, (2) 

the transition probability (|l]) transforms immediately into 

^{t,r^r') = W{t,ry)^{fJ') . (3) 

We remark that the transition rate W{t, r, r') depends only on the nature of the interaction 
between the test particle and the bath and doesn't depend on the population of the test 
particle in the starting and arrival sites. 

Kinetical interaction principle: In this paper, we will study kinetics coming from the 
transition probabilities @ and (P when the function 7 satisfies the condition 
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7(/, n ^ <f)_ 

where is a positive real function. This condition imphes that '-f{f,f')/K,{f) is a 

symmetric function. Then we can pose 7(/, /') = K{f)b{f)b{f')c{f, /') where b{f) and 
/') = c(/', /) are two real arbitrary functions. It will be convenient later on to intro- 
duce the real arbitrary function a{f) by means of 

«(/) ^ ^ . (5) 

and write 7(/, /') under the guise 

7{fJ') = aif)b{f)c{fJ') . (6) 

We claim at this point that 'y{f,f') given by (H) with a(/) and linked through 
is the most general function obeying the condition (|). We wish to note that 7(/, /') is 
given as a product of three factors. The first factor a{f) is an arbitrary function of the 
particle population of the starting site and satisfies the condition a(0) = because if the 
starting site is empty the transition probability is equal to zero. The second factor &(/') is 
an arbitrary function of the arrival site particle population. For this function we have the 
condition 6(0) = 1 which requires that the transition probability does not depend on the 
arrival site if, in it, particles are absent. The expression of the function b{f') plays a very 
important role in the particle kinetics, because stimulates or inhibites the transition r —>■ r', 
allowing in such a way to consider interactions originated from collective effects. Finally, 
the third factor c(/, /') takes into account that the populations of the two sites, namely / 
and /', can eventually affect the transition, collectively and symmetrically. 

The function 7(/, /') given by (^) defines a special interaction which involves, separately 
and/or together, the two particle bunches entertained in the starting and arrival sites. We 
observe that this interaction is different from the one depending on the coordinates of the 
sites involved in the transition which one takes into account by means of the functions T 
(cross section) in (|I]) or by W (transition rate) in (^. In order to explain the nature of the 
interaction introduced by the function 7(/, /') we start by considering the case 



7(/,/') = /(!-/') • (7) 



It is well known p0| , |36| that this particular expression for the 7(/, /') given by (J^) takes into 



account the Pauli exclusion principle and defines completely the fermion kinetics. Other 
expressions of the function 7(/, /') take into account interactions introduced by the gener- 
alized exclusion-inclusion principle |3^, the Haldane generalized exclusion principle p^ , |3T 



the Tsallis principle underlying the nonextensive statistics etc. We observe that the 

above mentioned principles impose the form of the collisional integral in the kinetic equa- 
tions through the choice of 7(/, /'). It is worth noting that in the cases of Haldane statistics 
the particular expression of 7(/, /') is originated from the fractal structure of the single par- 
ticle Hilbert space, its dimension depending on the particle number in the considered state 
||T6| , |3l|| . Also the Tsallis statistics is originated from the fractal structure of the relevant 
particle phase space p. 

Taking into account that particular choices of 7(/, /') reproduce the already known 
principles above mentioned, we can see the function 7(/, /') as describing a general principle 
which we call Kinetical Interaction Principle (KIP). The KIP defines a special collective 
interaction which could be very useful to describe the dynamics of many body systems. As 
we will see in the following sections, the KIP both governs the system evolving toward the 
equilibrium and imposes the stationary state of the system. 



III. KRAMERS GENERALIZED KINETICS 

In the following we study the particle kinetics in a 2n-dimensional phase space of a 
dilute system composed by identical particles interacting with an equihbrated bath. The 
procedure which we use in the present section to derive the evolution equation of the system, 
is a generalization to the non linear case of the standard procedure, involving the Kramers- 
Moyal expansion and the first neighbor approximation, which was introduced firstly to 
study the linear kinetics. We indicate with x and v the position and the velocity variables, 
respectively. The particles evolve under an external potential V = V{x). The evolution 
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equation for the distribution function f = fit, x, v) is given by 

^ = / [Tr{t,x,v' ^ v) - TT{t,x,v ^ v')] d^'v' 
at Jn 



(8) 



where df/dt is the total time derivative while the transition probability according to the 
KIP is given by 



nit,x,v ^ v') = Wit,x,v,v')^{f, f) . 



(9) 



Let us write the transition rate as W{t, x, v, v') = w{t, x, v, v' — v), where the last argument 
in w represents the change of the velocity during the transition. In the following, for sim- 
plicity, we indicate explicitly only the dependence on the velocity variables of the functions 
w{t, X, V, v' — v) and f{t, x, v). We start by writing Eq.(H) as follows 

df r 

-j-= / w{v + y,y) -f[f{v + y), f{v)]d''y 
dt Jn 

- I wiv,y)^[f{v),fiv-y)]d''y . (10) 



For physical systems evolving very slowly w{v,y) decreases very expeditiously as y 
increases and we can consider only the transitions for which v ± y ^ v. At this point we 
make use of the two following Taylor expansions 



w{v + y,y)-f[f{v + y), f{v)] 



OO -| 

E 



m=0 



d^{w{u,y)^[f{u)J{v)]} 



Z/oiZ/02 ■■■yam 5 



llfiv), fiv - y)] 



m=0 



ml 



g-7[/(^),/H] 
(9m„,(9m„2...9mo„ 



Z/ai2/a2 ■■■Vam 5 



and after substitution in (|T0|) we obtain the following Kramers-Moyal expansion. 



df{t,x,v) 
dt 



E 

m=l 



5'"{C«ia2...a„(i, X, U) 7[/(t, X, U), f{t, X, v)]} 



d^^[fit,x,v)Jit,x,u)] 

dUa^dUa2^^^dUam 



(11) 



where the m-th order momentum C«i«2-- am(^5 ^i of the transition rate is defined as: 

Caia2...a^{t,x,v) = ^ f ya,ya2- ■ -Va^w {t, X , V , y)d''y . (12) 
ml Jn 



We remark that from Eq. (|TT]), where the dependence on all the variables t, v is indicated 
explicitly, we can obtain as a particular case Eq. (7) of ref. ||3^ . 

In the frame of the first neighbor approximation only the first order (drift coefficient) Q 
and the second order (diffusion coefficient) C,ij momenta of the transition rate are considered. 
Indicating again explicitly only the dependence on the velocity variables Eq. (|Tl|) reduces to 
the following non linear second order partial differential equation 



dl 
dt 



9{0(n)7[/(n), /(^ ^ 0(^)- ^^""^^ 



dui 



duiduj 



dHf{v)J{u)] 



duiduj 



-' u=v 



which after taking into account the two identities: 



d{QiuUfiu),f{v)]} d^[fiv)Jiu)] 
H Ci[v)- 



d_ 



dui 



{avw{v)j{v)]}, 



dui 



and 



d^C.,{uUf{u),f{v)]} 



duiduj 



dHf{v)J{u)] 



duiduj 



-T^i 7[/(*^),/(^^)] 



+Cii(« 



dvj 

d^[fiu)Jiv)] d^[fiv)Jiu)] 



assumes the form 



dl 
dt 



dui 



d_ 

dvi 



du-i 



0(«) + 



dvj 



l[fiv)J{v)] 



d^[fiu)J{v)] d^[f{v)J{u)] 



duj 



duj 



(13) 



(14) 



Finally Eq. (|14D can be rewritten as 
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dt dvi 
with 7(/) = 7(/, /) and 



+ 



§^)7(/) + 0,7(/)A(/) 



di 



(15) 



A(/) 



5/ 7(/',/) 



-1/'=/ 



By taking into account the condition (|), the function A(/) simphfies as 

dlnKif) 



A(/) 



df 



(16) 



and Eq.(^) becomes 



(it dvi 



|7(/)+7(/) 



df 







(17) 



We assume the independence of motion among the n directions of the homogeneous and 
isotropic n-dimensional velocity space and pose Q = Ji, Qj = D6ij, being J and D the drift 
and diffusion coefficients, respectively. Moreover we introduce the function U by means of 



^dU _]_fj^dir 
dv D \ dv , 



(18) 



with j3 a constant. In the following we will consider the case where U = U{y) depends 
exclusively on the velocity. Taking into account that the potential V = V{x) depends only 
on the spatial variable, Eq. (|l^) can be written as 



^i^^^ = ^ (DivMf)-^[(3[Vix)+Uiv)-^^] + \nnif) 



with fi a constant and the expression of the total time derivative given by 



(19) 



d 



d ^ 1 dU{v) d 1 dV{x) d 



dt dt m dv dx m dx dv 

Equation (|19|) represents the evolution equation of the particle system in the Kramers picture 
and describe a non linear kinetics. This non linear evolution equation can be written in the 
form 



dt dv 



\. rr^d SIC 







(20) 
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where S]C/Sf is the functional derivatives of the functional /C defined through 



K = - I (TxcTv I dfln- 
Jn J K 



Us) 



where the stationary distribution fg = /(oo, x, v) is defined through: 



(21) 



In K,( 



-P[V{x)+Uiv)-fx\ . 



(22) 



We remark now that from Eq. (^) it follows immediately that the stationary distribution 
maximizes /C and can be obtained from a variational principle: 



(5/C 







df 







/ = /. 



Sf dt 
It is easy to verify that the functional /C increases in time 



die 



n of dt 

5K, d 
d'^xd^'v— — 
n of ov 



(23) 



In order to study the behaviour of the functional }C(t) when t — > oo we introduce the function 
(t(/) = — J df In K,{f ) so that K,{f) can be written as /«(/) = exp [—dcr/df]. Now we are able 
to calculate, in the limit t oo, the following difference 

/C(t)-/C(oo)= / d-xd'^v [a{f)-a{f,) + {f-f,)\nK{f,)] 
Jn 

Mfs) 



d'^xd'^v 



n 



d'^xd^'v 



n 



C^{f) - <fs) - {f - fs)- 



(24) 



.2 

and assume that d'^a{f) / df"^ < 0. This requirement is satisfied if the function n{f) obeys 
to the condition dn{f)/df > and consequently we have /C(t) < /C(oo), then /C assumes 
its maximum value for t = oo. The inequalities dK,{t)/dt > and /C(t) < /C(oo) imply 
that — /C is a Lyapunov functional and demonstrate the H-theorem. The functional /C is the 
constrained entropy of the system and results to be the sum of two terms: JC = S + Sc where 
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s 



J^rx^v J df In K{f) , (25) 



is the entropy of the system and Sc = —13{E — iiN). The energy E of the system is given by 

E = I d''xd''v[V{x) + U{v)]f . (26) 

We remark that, being k(/) an arbitrary function, the H-theorem has been verified in a 
unified way, for a very large class of non linear systems interacting with a bath. 

IV. BOLTZMANN GENERALIZED KINETICS 

In the diffusive approximation, adopted in the previous section to describe the changes 
of the particle states, the system is coupled with its environment. In this frame, meanwhile 
the particle system, interacting with the bath, evolves toward the equilibrium, both its 
entropy S and K, increases monotonicaly. In the present section we will consider the particle 
system isolated from its environment and describe its time evolution in the frame of the 
more rigorous Boltzmann picture. In presence of external forces derived from a potential 
the total time derivative is defined as 

d d d I dVix) d 
+ v- ^ ' 



dt dt dx m dx dv 
and the evolution equation assumes the form: 



dl 

dt Jn 



f d''v'd''Vid''v[[TT{t,X,v' V,V[ Vi) 

-TT{t,X,V^v',Vi^v[)] . (27) 



Eq. ( p7D describes a non linear generalized kinetics, the transition probabilities being defined, 
according to the KIP, as: 

7r(t, x,v ^ v', Vi v[) 

= Tit,x,v,v',v,,v[)^{fJ')^{f^J[) . (28) 
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In (|28|) we have posed f = f{t,x,v), f = f{t,x,v') and analogously fi = f{t,x,Vi), 
f[ = f{t,x,v[) in order to consider point-like binary collisions. A symmetry is im- 
posed to T = T{t, X, V, v', vi, v[) by the principle of detailed balance T{t, x, v, v', vi, v[) = 
T{t, X, v', V, v[, Vi), so that Eq. (pTf ) assumes the form 



dt Jn 



I T(t, X, V, v\ vi, v\] 

Jn 



x[7(/',/)7(/{,/i)-7(/,/')7(/i,/0] • (29) 
Alternatively, by taking into account one can write: 

f=fd'^v'd^v,d'^v[Tc{f,f')c{fiJi) 
at Jn 

X [a{f)b{f)a{f[)b{f,) - a{f)b{f)a{f,)b{f[)] . (30) 
We note that when c(/, /') = c(/i, f[) = 1, Eq. (0) reduces to the equation recently con- 



sidered in ref. |3J]. In the following we will describe the system using Eq. (^91) which can 
be rewritten as 



dt Jn 



f d-v'd^v,d-v[ T Q 7(/', /) 7(/{, /i) , (31) 
Jn 

where the auxiliary function Q = Q{f, f, fi, f[) < 1 is defined as 



^ lifJ'hifiJD 
lif'JhifiJi) ■ 

Taking into account the condition the function Q becomes Q = 1 

and can be written immediately in the following form 

Q = 1 - exp [ In Kif) + In - In - In «:(/()] • 



We consider now the system at the equilibrium. From the evolution equation (^) we have 
df /dt = 0, and Q = 0. The first condition df /dt = taking into account the definition of 
the total time derivative, implies for the stationary distribution f = fs that fs = /s[mi7^/2 + 
y(a3)]. On the other hand from the condition Q = we have 

In + In - In «:(/;) - In K{f:,) = . 
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This last equation allows us to conclude that the quantity ln/t(/s) is a collisional invariant 
for the particle system. If we suppose that the binary interparticle collisions conserve the 
particle number and the kinetic energy 

1 2,1 2 1 /2,1 /2 

-mv H — mVi = -mv H — mv^ , 

we have that the quantity — /5 ^mv^ + V{x) — is the more general collisional invariant. 
Then we obtain the condition 

\nK{f,) = -l3^^mv^ + V{x)-fi^ , (32) 

which defines the stationary distribution, so that Q becomes 

Q= 1 — exp 



4f.) 4f,i) 4fl) 

After introducing the functional /C by means of (pTD with fs is given by (|32D the quantity Q 
can be written in the form: 



5f 5h 5f' 6f[ 
Finally the evolution equation (|3lD assumes the form 



dl__ 
dt Jn 



(Tv'd-v.d-v'.T 7(/',/)7(/(,/i) 



(33) 



5f 5h 5f 6f{ )] 
From the structure of (p3|) we have that the stationary distribution fs = f{oo,x,v) can be 
obtained from a variational principle 



sk: 







df 



; / = /. . 



Sf dt 

We study now how the functional /C evolves in time. Its time derivative is given by 



die 



<f) df 



d^xd^v lu- 
ll K[fs) dt 

I d''xd''vd''v'd''vid''v[ T In 
Jn 



m ' f^ifs) 

x[7(/',/) 7(/{,/i)-7(/,f) 7(/i,/0] 



(34) 
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The symmetry of the integrand function permit us to write (|3^) as it follows 

. 1 ^ 



die 

Hi' 



- 1 <rx<rvrv'<rvi<rv[ ^ t 



In- 



<fsi) «:(/;) 
X [7(/',/)7(/{,/i)-7(/,/')7(/i,/0] • 

After taking into account the expressions of Q we have: 

, 1 



(35) 



die 

It' 



I d''xd''vd^v'd'^vid''v[ - T 
Jn 4 

[-Qln(l-Q)] 7(/',/) l{f[Ji] 



(36) 



We observe now that Q < 1 and then we have — Qln(l — Q) > 0. This implies that the 
integrand function in (^) is a non negative function. Then we conclude that dJC/dt > 0. 
Starting from the definition of /C and following the procedure adopted in Sect. Ill we obtain 
/C(t) < /C(oo). We can conclude at this point that — /C is a Lyapunov functional. It is easy 
to verify that 



JC = S-f3{E- fiN) 



(37) 



where the entropy of the system 5* is defined through (^) while its energy E is given by 

1 



E= I d'^xd^'v 
In 



-mv + V{x) 



f , 



(38) 



which is a conserved quantity, dE/dt = 0, as the particle number A^. Consequently ( p7D can 
be written as 



]C{t) = S(t) + constant. 

The H-theorem for the isolated nonlinear system follows immediately: 

dS 



(39) 



dt 



> 



S{t) < S{oo) 



(40) 
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V. SOME KNOWN STATISTICS 



In this section we will show that the formalism previously developed permit us to con- 
sider, in a unitary way, the already known statistical distributions. For simplicity we will 
discuss the case of distributions depending exclusively on the velocity. Firstly we observe 
that the stationary distribution /, defined through «;(/) = exp(— e) with e = (5{mv'^ /2 — /i), 
can be obtained as steady state of a Fokker-Planck (FP) equation describing the kinetics of 
brownian particles for which results U — mv^ /2 and D — const. The same distribution can 
be viewed as steady state of a Boltzmann equation, describing free particles interacting by 
means of binary collisions, conserving the particle number, momentum and energy. In this 
section we will write the evolution equations (FP and/or Boltzmann) of some distributions 
available in literature to illustrate the relevance of the approach adopted in the previous 
sections describing the non linear particle kinetics. 

Maxwell- Boltzmann statistics: We start by considering the MB statistics given by / = 
Z~^exp(— e). It is readily seen that the related kinetics is defined starting from a(/) = /, 
h{f') = 1, while the symmetric function c(/, /') remains arbitrary. Then we have «(/) = / 
and 7(/) = /c(/). In the Boltzmann picture the evolution equation becomes 



In the simplest case c(/, /') = 1 we obtain the standard linear Boltzmann and FP equations. 
We observe that there is an infinity of ways (one for any choice of c(/, /')) to obtain the MB 
distribution. 

Bosonic and fermionic statistics: We consider now the case of quantum statistics namely 
the Fermi-Dirac (77 = —1) and Bose-Einstein (77 = 1) statistics defined by means oi f — 
Z~^(expe — r])~^. The kinetics now is defined through a{f) = f and &(/') = 1 + rjf while 




(41) 



while, in the Fokker-Planck picture we have 




(42) 
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again the function c(/, /') remains arbitrary. We have consequently = //(I + r]f). In 
the Boltzmann picture the evolution equation becomes 
Of 



= |^rfVrf>dX Tc(/,/')c(/i,/0 

X [/'(I + vf)fia + vfi) - /(I + + rim 



(43) 



and reduces to the well known Uehling-Uhlenbeck equation if we choose c(/, /') = 1. This 
choice for c(/, /'), in the frame of the FP picture, implies 7(/) = /(I + rjf) and we obtain 
the following evolution equation PB|| : 



dt dv 



D(3mvf{l + r]f) + D 



dv 



(44) 

In order to show that a kinetics, different from the one described by (PI), also reproducing the 
bosonic and fermionic statistics, exists, we consider c(/, /') = [1 + 1]^/ f f')~^ or alternatively 
c{f,f') = [1 + vif + It is easy to verify that in both the cases we have c(/) = 

(1 + i]f)~^ and 7(/) = /. Now the evolution equation in the FP picture becomes ||1 1|| : 

df 



dt 



d_ 
dv 



1 d 

D(3mvf + D- — \Yiil + r]f) 
7] dv 



(45) 



Intermediate statistics: The quantum statistics interpolating between the bosonic and 
fermionic statistics has captured the attention of many researchers in the last few years. 
A first example of intermediate statistics can be realized by considering in the distribution 
/ = Z~^(exp e — r/)"^, previously examined, the parameter rj as being continuous: < r/ < 1. 
For ?7 7^ ±1 we have a quantum statistics different from the Bose or Fermi statistics. A second 
intermediate statistics is the boson-like (+) or fermion-like (— ) quon statistics which can 
be obtained easily by posing a(/) = [/] ^ and &(/') = [1 ± /'] where [x] ^ = {q^ — g~^')/2 In q 
and g G R. If we choose for simplicity c(/) = Cg = 21ng/(g — q^^), the evolution equation 
in the FP picture becomes |p2|| : 



A third intermediate statistics is the Haldane-Wu exclusion statistics which can be obtained 
starting from the kinetics defined by setting a(/) = / and &(/') = {l — gf Y [1 + (1 — (7)/'] 



(46) 



with < ^ < 1 |3T 
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Tsallis statistics: We consider the non extensive termostatistics introduced by Tsallis 
0. The relevant distribution / = Z^^[l — (1 — q)eY^^^~'^^ can be obtained naturally starting 
from the kinetics defined through lnfi;(/) = (/^^"^ — 1)/(1 — g) = In, /. The Boltzmann 
equation (|30| ) becomes now 



df 



dt Jn 



d-v'd-v^d-v[ T c{f, f)c{f,, f[)h{mh)h{f)h{f[) 



X 



exp (In, /' + In, f[) - exp (In, / + In, /i 



(47) 



In the FP picture the evolution equation is given by 



dl^d_ 
dt dv 



D^{f) (3mv + f 



dv 



(48) 



which for 7(/) = / reduces to the one proposed in ref. [pi 



VI. THE fv-DEFORMED ANALYSIS 

In the previous section we have considered some statistical distributions, quantum or 
classical, already known in the literature, depending on one continuous parameter. We will 
now turn our attention to the distribution / = Z~^(expe — r^)^^. We note that this quantum 
distribution can be viewed as a deformation of the MB one, which can be recovered as the 
deformation parameter t] approaches to zero. Another classical distribution which can be 
obtained by deforming the MB one, is the Tsallis distribution / = Z^^[l — (1 — g)e]^/(^^''\ 
The MB distribution emerges again as the deformation parameter g — >■ 1. 

In the present section we will study the main mathematical properties of a new, one 
parameter, deformed exponential function, while in the next section we will consider the 
induced deformed statistics. The deformed exponential is indicated by exp^^^(x), where k 
denote the deformation parameter, and we postulate that it obeys the following condition: 

exP(^j(x)expj^j(-a;) = 1 . (49) 

We start by observing that any function A{x) can be written in the form A{x) = Ae{x) + 
Ao{x) where Ae{x) = Ae{—x) is an even function and Ao{x) = —Ao{—x) an odd one. The 
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condition A{x)A{—x) = 1 allows us to express Ae{x) in terms of Ao{x) by means of Ae{x) 



1 + Ao{xy and consequently write the function A{x) in the form: A{x) = yl + Ao(x)^ + 
Ao{x). At this point it is obvious that a deformed exponential exp^^j(x) obeying (^9[) and 
depending only on one deformation parameter k, so that exp^^^ (x) ~^ exp x, can be written 
as 



exp^^j(x) 



1 



(50) 



or alternatively as 



exp^^^(x) = exp (^iarcsinh5(«(x)j . (51) 



In (|SD|) the generator (x) of the deformed exponential is an arbitrary function depending 
on the parameter k and obeying the conditions: 

Qni-x) = -5'k(x) ; g^{x) ~ kx . (52) 

Since it will be useful later on, we introduce the inverse function of exp^^^ (x) indicated by 
ln^^^(a;) and defined through exp^^^ In^^j(x) = x. It is easy to verify that 



In^^^ix) = g-' j , (53) 

where g~^{x) is the inverse function of g^ix). We remark that the deformed exponential 
can be defined by fixing the expression of the generator g^ix). We note that by choosing 
gK.{x) = smhux we can generate the standard undeformed exponential exp^^^(a;) = expx 
and logarithm In^^j (x) = In x functions. 

The K- exponential: In the following we consider the simplest deformed exponential (in 
following called ^-exponential), which is generated from g^ix) = nx and is given by 

1/k 



exp^^j(x) = [yl + K^x^ + KXj , (54) 
or equivalently by 

exp^^^ (x) = exp arcsinh kx] . (55) 



K 
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Obviously, we have exp^^^(x) = expx and for Vx e R it results exp^^^(x) e R"*". Furthermore 
we have exp^^j(O) = 1 and exp^_^^(a;) =exp^^j(a;), so that we can consider simply that 
K e R"*". A relevant property of exp^^^(a;) is that for Va e R 

Concerning the asymptotic behaviour of ^-exponential we easily obtain that 
exp (x) ~ 1 2ra I and exp (x) ~ |2/txh^/''''. We observe that the deformed ex- 
ponential is an increasing function dexp^^^{x)/dx > 0, e R. Finally, its concavity is 
(P exp^^^(x)/dx^ > for |«:| < 1 while for \k\ > 1 we obtain exp^^^{x)/dx'^ > ior x < Xc 
and d"^ exp^^^(x)/dx^ < ior x > Xc being Xc = (k^ — k,^)~^^^. 

The K-sum: Of course if we take into account that arcsinh x + arcsinh y 
— arcsinh {x\/l + + y\^T+~x^) we arrive immediately at the following relationship: 

exp^^j (x) exp^^^ (y) = exp^^^ © y) , (57) 

where the K-deformed sum or simply re-sum x ^y, is defined through 



X (B y = xy 1 + K'^y'^ + yvl + Kp^ . (58) 

We note that the re-sum obeys the associative law {xS)y)®z = x(B{y(Bz), admits a neutral 
element x®0 — 0®x — x and for any x exists its opposite x ® {—x) — 0. Moreover the 
commutativity property x ^ y — y ^ x holds, so that the real numbers constitute an abelian 
group with respect the re-sum. Since it will be useful later on, we define the re-difference as: 
xQy = X® {-y). 

The K-trigonometry: Firstly we define the re-deformed hyperbolic sine and cosine, starting 
from exp^^j(±a;) = cosh^^j(a;) ± sinh^^j(a;) which is the re-Euler formula and obtain 



2 
1 

cosh, , (x) — - 



-x\ 



exp^^j(x) + exp^^j(-x) 



It is straightforward to introduce the re-hypcrbolic trigonometry which reduces to the ordi- 
nary hyperbolic trigonometry as re ^ 0. For instance, the formulas: 
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cosh^ , (x) — sinh^ (x) — 1 , 

{K} ^ ' {K} ^ ' ' 

sinh^^j {x © y) +smh^^^ {xey) = 2 sinh^^^ (x) cosh^^j (x) , 

+ ^ r I. r \ sinh {x®y) 
tann, , (x) +tann. , I'y) = — ; ^ -— , 

and so on, still hold true. All the formulas of the ordinary hyperbolic trigonometry still hold 
true after expediently deformed. The deformation of a given trigonometric formula can be 
obtained starting from the corresponding undeformed formula, making in the argument of 
the trigonometric functions the substitution x + y ^ x ^y and obviously nx ^ x ^ x... ^ x 
(n times). 

The «;-De Moivre formula involving trigonometric functions having arguments of the type 
rx with r e R, assumes the form 

[cosh^^j(x) ±sinh^^^(x)]''= cosh^^/^j(ra;) ±sinh^^/^j(rx) . 

The K-derivative: It is important to emphasize that we can naturally introduce a re- 
calculus which reduces to the usual one as the deformation parameter re — > 0. In fact we 
can define the re-derivative of course through 

where the re-differential d x^^^ = limy^^. x Q y, after taking into account the identity 
{x ® y){x Q y) — x"^ — y^, is given by 

7 1. ^■ ~ dx 

d x,^, = lim X Qy = lim ^ = , = . 

The function x^^y can be easily calculated after integration 

x^^j = — In (^vT+~re%^ + Kx^ = — arcsinh rex , (60) 

and satisfies the relationship x^^y ~ i^®y){K}- observe that it is possible to write 

the re-derivative in the form 
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^ = ViT^^-^, (61) 

from which it appears clearly that the /^-calculus is governed by the same rules of the ordinary 

one. For instance, we can write 

dexp (x) dsinh (x) 

= exp^ (x) : ; — — = cosh, (x), 

and so on. We observe now that the introduction of the function x^^^ permit us to write 
exp^^^(x) = exp(^x^^y^ and analogously for the K-deformed trigonometric functions. For 
instance, we have sinh^^j(x) = sinh(^x^^j^, coshj^j(x) = cosh (^x^^j^ etc. This property of 
the K-deformed exponential and trigonometric functions permits us to consider their Taylor 
expansion in terms of the function x^^^ . For instance, it holds 

^^P^J^)=E^{^,.T ■ (62) 

m=0 

The K-integral: Starting from the definition of ^-derivative, we define the K-integral 
through 

f(x)dx,^, = I -jIM=dx . (63) 

All the standard rules of the undeformed integral calculus still hold true. 

The K- cyclic functions: The K-deformed sine and cosine can be defined through 
sin^^j(x) = — i sinh^^j(zx) and cos^^j(x) = cosh^^j(zx), respectively. It results: sin^^j(x) = 
sin(x,. ,) and cos, , (x) = cos(x,. ,) being x,. , = k"^ arcsin(Kx). Of course the ^-cyclic 
trigonometry, naturally can be introduced. 

The K-logarithm: We study now the inverse function of the ^-exponential namely the 
ft-logarithm which, after remembering its definition (|53|) and the expression of the generator 
gn{x) = Kx, can be written as 

From ( |64D we can see that Vx G =^ In^^j(x) G R and hi^^^{x) = Inx. Furthermore we 
have Inj^j(l) = and ln^_^j(x) = In^^j(x). Two relevant properties of the K-logarithm are 
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\n^^^{xy) = Inj^j(x) © \n^^-^{y) and ln^^j(x™) = mln^^^j(x), while its asymptotic behaviour 



is In (x) ~^ 



a;— >H-oo 



-oo. 



2k\ I**' — s> — oo and In, , ix) 

I {K}\ J 

In closing this section we consider the relationship linking the functions exp^^^(x) 
and Inj^j(x), here introduced with the already known in the literature Tsallis expo- 
nential expg(x) = [1 + (1 — g)x]^/*^^~'^^ and its inverse function, the Tsallis logarithm 
x^~'^ — 1)/(1 — g) namely: 



lng(x 



exp^^j^xj 



exp 



X + 



Vl + /t2x2 - 1~ 



In, , (x) + In, (x) 



VII. THE K-DEFORMED STATISTICS 

In this section we will consider a new statistical distribution, just as a working example 
of the theory presented in the previous sections. We start by considering a particle system 
in the velocity space and postulating the following density of entropy: 

^M) = -Jdf ln,^,{af) , (65) 

where the real constant a > is not specified for the moment. Eq. (^), for k ^ 0, 
gives the standard Boltzmann-Gibbs-Shannon density of entropy if we pose a = 1. We 
note that the above definition of ct^(/) implies that k(/) is related to the K-logarithm, 
through lnK(/) = \n^^^{af). It is immediate to see that for V/t G R, d?'a^{f)/dp < 0, 
independently on the value of a. Then for the entropic functional /C, we have dlC{t)/dt > 
and /C(t) < /C(oo). The entropy of the system given by S^=j^d'^va^{f) assumes now the 
form 

S, = -— I d% f ^ _ ^ f^A , (66) 

and reduces to the standard Boltzmann-Gibbs-Shannon entropy = — Jud'^v [ln(a/)— 1]/ 
as the deformation parameter approaches zero. This K-entropy is linked to the Tsallis entropy 

(T) 

5*^ through the following relationship 
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K 



1 a'^ (T) 1 a (T) 

2 l + K 1+^^ 2 1-K I-'' 



+ const. 



(67) 



Fzrst choice of a: We discuss now the case a = 1. The stationary statistical distribution 
corresponding to the entropy can be obtained by maximizing the functional /C 



which gives the standard classical statistical distribution as k — > 0. 

Second choice of a: Before introducing the second choice of the constant a we consider 
briefly the concept of entropy. It is to be understood that the entropy is an absolute measure 
of information only in the case of an isolated system, where the particle number and energy 
are conserved. In the case of a system interacting with a bath, only the particle number 
is conserved and the free entropy can't be used as an absolute measure of information. 
For this reason an entropy constrained by the particle number conservation and by the 
relevant energy mean value must be introduced. The constrain introduced by the particle 
number conservation is a special one, and because of its presence both in the case of free and 
interacting systems, can be taken into account directly in the definition of the entropy. For 
instance, for the linear kinetics, the stationary distribution / = exp(— /?[/) with partition 
function given by Z = Jud^'v exp(— /5f/), can be obtained using the variational principle, 
namely 6[Sq— P J-j^ d^v Uf] =0, where the functional 5*^ = —Jj^ d^v [ ln(Z/) — 1] / depending 
on the constant Z is the above mentioned constrained entropy, the particle number (in 
following we pose = 1) being conserved. 

It is clear that, for analogy, also in the case of the non linear kinetics with k 7^ 0, we 
can choose a = Z in the expression of S*^ given by (|^), so that the stationary statistical 
distribution 




(68) 



One arrives to the following distribution 




(69) 




(70) 
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with Z = J'jid'^v exp^^j (— /9?7) can be obtained by considering the following variational 
principle 



6 



S^-P I (TvUf 
n 



(71) 



Of course, the expression of / depends on the potential U and, in the particularly inter- 
esting case of Brownian particles, U = mv^/2, after tedious but straightforward calculations 
we can write the distribution function as: 



f- 



Pm\K\ 



n 



IH — nlK] 



, s exp, I mv 



(72) 



where n = 1, 2, 3 is the dimension of the velocity space and \k\ < 2/n. The distribution given 
by ([7^) reduces to the standard Maxwell-Boltzmann one / = (/?m/27r)"'/^ exp(— /3mf ^/2) as 
the deformation parameter k approaches to zero. This limit can be performed easily by 
taking into account the Stirling's formula: r(2) ~ -\/27r exp(— z), holding for 

We write now the evolution equation, whose stationary state is described by ([721). After 
indicating with f{t,v) the time depending statistical distribution, which for t— >oo reduces 
to (|72D , we introduce the new function p = Z^of, and remember that a = Z^. The evolution 
equation of the function p, in the Fokker-Planck picture, in the case of Brownian particles, 
and by choosing for simplicity 7(/) = /, assumes the form 



(9p 
'di 



d_ 
dv 



/ K -K^ 9p 



(73) 



where k = Dm(3oo- In the Boltzmann picture the evolution equation is structurally similar 
with the one of Tsallis kinetics. The only difference is that the Tsallis logarithm is replaced 
now by the K-logarithm. 

We conclude this section by considering a new quantum distribution, describing particles 
with a behaviour intermediate between bosons and fermions, which can be constructed 
starting from ([691). We impose for the entropy density the following expression: 



(74) 
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being t] a real parameter. After maximization of the constrained entropy or, equivalently, 
after obtaining the stationary solution of the proper evolution equation associated to (fTif), 
one arrives to the following distribution 

/ = — • 

exp^^j(e)-r/ 

We note that (|75D for rj = reduces to (^) while for 77 = 1 becomes a K-deformed Bose- 
Einstein distribution and for 77 = —1 becomes a K-deformed Fermi-Dirac distribution. Fi- 
nally for 77 7^ 0, ±1 the ( [75|) becomes a quantum distribution which defines an intermediate 



statistics and can describe anyons likewise of the distribution considered in ref. [Tgj3^. Ob- 
viously, other quantum intermediate statistical distributions (generalized Haldane statistics, 
generalized quon statistics, etc.), can be constructed starting from (|69D. 

VIII. APPLICATIONS OF k-DEFORMED STATISTICS 

It is worth to note here that k— deformed statistical distribution given by ( [69| ) and ( |70D 
is obtained by extremization, under constraints, of the entropy S*^ which is given by ( |67D 
in terms of Tsallis entropy. This interesting result permits us to use k— deformed statistical 
distribution to study physical systems where we can find experimental evidences for the 
physical relevance of the Tsallis entropy like, for instance, in 2d turbolent pure-electron 



plasma |p^, solar neutrinos |38-40|, bremsstrahlung |41], anomalous diffusion of correlated 



and Levy type p7| , |42[| , self-gravitating systems cosmology among many others. In 



the following we consider two examples of application of /t-deformed distribution, just to see 
the values the parameter k assumes. 

In ref. [^^-^ , firstly, the problem of the solar neutrinos is considered and a solution, in 



the frame of the non extensive statistics is proposed. It is well known that the solar core, 
where the energy is produced, is a weakly non ideal plasma. In fact, density and temperature 
condition suggest that the microscopic diffusion of ions is non standard: diffusion and friction 
coefficients are energy dependent, collisions are not two-body processes and retain memory 
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beyond the single scattering event. For this reason the equihbrium velocity distribution of 
the ions is slightly different in the tail from the Maxwellian one, as argued also by Clayton 



45|,|46[|. Consequently the reaction rates are sensibly modified and, at last, the neutrino 



fluxes which are experimentally detected, have calculated values different from the standard 
ones. With the hypothesis that the velocity distribution of the ions in the solar core is a 



K-deformed distribution, we can reproduce the experimental data, analyzed in ref. ||38|-^ 
when K = 0.15. 



In ref. it is shown that the observational data concerning the velocity distribution 
of clusters of galaxies can be fitted by a non extensive statistical distribution. If we adopt 



the K-deformed statistical distribution to analyze the data reported in ref. we obtain a 
remarkable good fitting when n = 0.51. 



IX. CONCLUSIONS 

In the present effort we have studied from a general prospective the kinetics of non linear 
systems in the frame of two pictures, namely Kramers and Boltzmann. The results, can be 
summarized as follows: 

The KIP governs the time evolution of the non linear system by means of the func- 
tion 7(/, /') = a(/) b{f') c{f, /') and imposes its steady state through the function = 
a{f)/b{f). The steady state fs can be obtained as stationary solution of its evolution equa- 
tion. 

The KIP imposes the entropy form of the non linear system, which is given by ( ^5|) both 
in Kramers and in Boltzmann pictures. The constrained entropy /C, given by (0), satisfies 
the H theorem when dK,{f)/df > and the fs can be obtained also from the maximum 
entropy principle for /C. 

In the case of Brownian particles where U{v) = |mf^, the stationary state fs and then 
the constrained entropy /C assume the same values both in Kramers and in Boltzmann 
pictures. Being K{f) an arbitrary function, the H-theorem has been demonstrated in a 
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unified way for a very large class of isolated or interacting with a bath non linear systems. 

Finally, we have considered within the formalism here developed, some statistical distri- 
butions already known in the literature. On the other hand, as a working example of the 
theory here presented, we have introduced the new K-deformed statistics. After studying 
the main properties of the new statistics, we have discussed two applications to real physical 
systems. 
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